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Abstract. We used red clump stars to measure the surface mass density of the Galactic disk in the solar neigh- 
bourhood. High resolution spectra of red clump stars towards the NOP have been obtained with the ELODIE 
spectrograph at OHP for Tycho-2 selected stars, and nearby Hipparcos counterparts were also observed. We de- 
termined their distances, velocities, and metallicities to measure the gravitational force law perpendicular to the 
Galactic plane. As in most previous studies, we applied one-parameter models of the vertical gravitational po- 
tential. We obtained a disk surface mass density within 1.1 kpc of the Galactic plane, Si.ikpc — 64 ± 5 Mqpc -2 , 
with an excellent formal accuracy, however we found that such one-parameter models can underestimate the real 
uncertainties. Applying two-parameter models, we derived more realistic estimates of the total surface mass den- 
' sity within 800 pc from the Galactic plane, Eo.8k P c = 57-66 Mqpc -2 , and within 1.1 kpc, Ei.ikpc — 57-79 Mqpc -2 . 

This can be compared to literature estimates of ~40 Mqpc -2 in stars and to 13 Mqpc -2 in the less accurately 
measured ISM contribution. We conclude that there is no evidence of large amounts of dark matter in the disk 
and, furthermore, that the dark matter halo is round or not vey much flattened. 

A by-product of this study is the determination of the half period of oscillation by the Sun through the Galactic 
plane, 42 ± 2Myr, which cannot be related to the possible period of large terrestrial impact craters ~ 33-37 Myr. 

Key words. Stars: kinematics - Galaxy: disk - Galaxy: fundamental parameters - Galaxy: kinematics and dynamics 
- Galaxy: structure - solar neighbourhood 
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1. Introduction thereby allowing one to constrain the determination of 

■ . . . the total volume mass density in the solar neighbourhood, 

I his paper is the extension of the previous works pub- n , „ , _ , , . , 

, — — — , . „ tn i now called the Oort limit, iwo decades ago, a seminal 
hshed by Soubiran et al. 2003, hereafter Paper 1 and , ,. n , „ rrrom i u 

Siebert et al (2003 hereafter Paoer II) which orobed the lm P rovement was achieved by Bahcall 0223, who built 
oieueii ei ai. zuuo iiLicaiiei _rapei ill, wiiicn pioueci Hie . , . ^ , . .. n n n n . . 

' — , , r , r, a consistent Galactic vertical potential linked to current 

properties of red clump stars withm 100 pc of the Sun . , , , . . , , , . 

, ,. , „ , . „ , knowledge of the kinematics and the density distributions 

and at larger distances towards the North Galactic Pole c . , ,. „. , , , rmrrm c n j 

„ mT1 , TTr ° ... , r , , of stellar populations. Bienayme et al. (1987) followed a 

NGP . We obtained a new determination of the local sur- . ., , , „ , 

\ ' , . ,„ .. . . very similar approach by constraining the Galactic verti- 

face mass density. We discuss this new result with respect , . , , , , , , , ^ , . , , , 

, J . . , cal potential through global Galactic star counts and the 

to recent works, and we comment on the different ways the , , , , ,, , ,, , , . , . 

current knowledge on the stellar population kinematics. 

uncertainties have been estimated and the resulting con- 
sequences for the estimated mass density of the Galactic Later, a ma jor step forward was made by Kuijken & 

disk- Gilmore l|1989[) with a new sample of K dwarfs towards 

SinCe th l P i2 I i!!^ g WOrks ° f Ka P te y n tUEia and, thc South Galactic Pole (SGP) tracing the vertical poten- 

later, Oort H1932J 11960J) , regular improvements have been tiaL They used the same gtars to measure both the ver tical 

obtained in determinin g thc vertical Galactic potential, density distribution and kinematics. This had thc immc- 

* Based on observations made at the Observatoire de Haute dia te consequence of considerably reducing uncertainties 

Provence (France). Based on data from the Hipparcos astrom- existing in previous works, where different samples were 

etry satellite. used to determine both the vertical density and kinemat- 
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Thereafter, regular advances occured with improved 
stellar samples and accurate corrections of systematic ef- 
fects by Flynn and collaborators ljT3§lllT3$8l 12UUU1 !2UUl|l . 
in our Papers I and II, and by other authors (for instance 
Korchagin et al., 2003). 

A decisive moment was, of course, the arrival of 
Hipparcos observations ( |ESA, 1997| ) that allows precise 
calibration 1 of stellar parallaxes and absolute magnitudes. 
This has allowed robust estimation of distances and also 
useful measurement of tangential velocities. An immedi- 
ate application has consisted in probing the potential close 
to the Galactic plane within 100-200 pc directly. Likewise, 
Hipparcos data gave immediate access to the Oort limit 
(PhamUnnHI Creze et al. H998al[l998bl Holmberg & Flynn 
Korchagin et al. lMiSjl . 

In this general context, our paper describes the obser- 
vational extension of the red clump samples analysed in 
Paper II, and gives a new dynamical determination of the 
total local surface mass density. We also discuss what are 
probably the real current constraints obtained on the sur- 
face mass density, and comment on other results obtained 
in previous papers. 

There is no perfect agreement yet between the var- 
ious recent determinations of the vertical potential per- 
pendicular to the Galactic plane. Samples remain small; 
methods and analysis are probably not yet optimized; and 
some assumptions, like full phase mixing and stationarity, 
are difficult to check. Furthermore, useful and complemen- 
tary information are not used optimally like the change of 
metallicities with kinematics. 

Even if many systematic effects can now be conve- 
niently considered and corrected, the lack of large un- 
biased stellar samples with radial velocities prevent us 
from examining the stationarity of stellar tracers in de- 
tail, which is a central question that will need further 
examination. We may, however, note that the Hipparcos 
proper motions and tangential velocities have allowed the 
3D velocity distribution in the solar neighbourhood to 
be probed. The phase mixing appears to be 'slow' for 
horizontal motions within the Galactic plane. The corre- 
sponding period of the epicyclic motions is 169 Myr, and 
a few streams are still visible in the (w, v) velocity space 
(Chereul et al. 115991 and Famaey et al. !2UUl|) . For vertical 
motions, the oscillation period is shorter, 86 Myr, or half 
the epicyclic period, and only one velocity stream is still 
clearly visible (associated with the Hyades cluster); other- 
wise the (z,w) phase space corresponding to the vertical 
motions seems to be phase-mixed. 

Future advances are expected with the measurement 
of the disk surface mass density towards regions away 

1 We note that the correction of systematic effects should be 
easier with Hipparcos data, since the distribution of errors is 
understood so well. However, the way the correction of system- 
atic effects is applied always depends on the astrophysical ques- 
tion examined. The bias of Lutz-Kelker, Malmquist or others 
must be cautiously considered to achieve a proper correction; 
see for instance a discussion by Arenou et al. £1999 1 . 



from the solar neighbourhood. Local kinematics still car- 
ries non-local information about the structure of the disk: 
for instance, the coupling between the vertical and (u, v) 
horizontal motions of stars in the solar neighbourhood is 
directly linked to the scale length of the total mass distri- 
bution in the disk (Bienayme, 1999). 

Analysis of stellar IR surveys, 2MASS, DENIS or 
DIRBE (Smith et al.. l2004|) will allow minimization the ef- 
fects of the extinction and uncertainties on distance scales. 
Surveys like the RAVE project (Steinmetz, 2003) will in- 
crease the number of available radial velocities by one or 
two orders of magnitude: a few hundred thousand bright 
stars with I < 12. The next step, an ESA cornerstone 
mission, will be the GAIA project (Perryman et al.. 12001). 
but new methods of analysis should be prepared. Classical 
analyses, like the one applied in this paper, would certainly 
be insufficient for fully investigating the huge amount of 
data expected. 

In Sect. 2 we describe selection of the 3 samples that 
we use: a local one and 2 distant ones towards the North 
Galactic Pole, optimized to include a large fraction of 
clump giants. Section 3 is devoted to methods and ex- 
plains how we determinated of the vertical potential and 
the disk surface mass density. Our discussion and con- 
clusions are given in Sects. 4 and 5. In Soubiran et al. 
(2005, hereafter Paper IV), we describe the improvement 
of our local and distant samples in detail as compared to 
Papers I and II, and we analyse these samples in terms of 
the properties of thin and thick disk populations. 

2. The survey 

To determine the vertical force perpendicular to the 
Galactic plane, we measured the spatial and the verti- 
cal velocity distributions of a test stellar population as a 
function of vertical height. As far as possible, this test pop- 
ulation must be homogeneous and unbiased with selection 
criteria that are independent of velocity and distance. It 
must also be in a stationary state. For this purpose, we 
used one local and two distant samples of red clump gi- 
ants selected within the same B — V colour window and 
the same absolute magnitude window; our selected NGP 
clump stars are the distant counterparts of our selected 
Hipparcos stars. We find that at magnitude V — 9.5 to- 
wards the NGP, half of stars, with 0.9 < B - V < 1.1 
are clump stars while at magnitudes V < 7.5, more than 
80 percent are clump stars. Redder and bluer clump stars 
do exist. The blue cut removes the most metal poor stars 
([Fe/H] < —0.60) very efficiently but allows us to reach 
fainter magnitudes with low contamination by main se- 
quence stars. The red cut was applied to minimize the 
contribution of subgiant stars and of other giants on their 
first ascent of the giant branch. 

The distant sample is the extension to larger distances 
from the Galactic plane of the NGP sample that was pre- 
viously analysed in Papers I and II. It was built from a 
preliminary list of red clump candidates from the Tycho-2 
star catalogue (H0g et al. 120001) . High resolution spectro- 
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scopic observations were used to confirm the red clump 
stars, to separate them from the other stars, and to mea- 
sure radial velocities. We also improved the local sample 
of 203 red giants by measuring new radial velocities and 
metallicities for 88 of these stars. The selection, observa- 
tion, and reduction of the two samples is briefly described 
below and explained in Paper IV. 

2.1. The Hipparcos red clump stars 

The local sample of 203 nearby red clump giants was se- 
lected from the Hipparcos catalogue according to the fol- 
lowing criteria : 

7r > 10 mas 
<5icrs > -20deg 
0.9 < B-V < 1.1 
< My < 1.3 

where ir is the Hipparcos parallax and <5icrs 2 the decli- 
nation. The Johnson B — V colour was transformed from 
the Tycho-2 B^ — V? colour by applying Eq. 1.3.20 from 
|ESA, 1997| 

B-V = 0.850 (B T -V T ). 

The absolute magnitude My was computed with the 
V magnitude resulting from the transformation of the 
Hipparcos magnitude H p to the Johnson system with the 
equation calibrated by Harmanec (1998). 

We searched for radial velocities and metallicities for 
these stars in the literature. Our source of spectroscopic 
metallicities is the [Fe/H] catalogue (Cayrel et al.. l2001|) . 
in which we found a fraction of our stars in common with 
McWilliam ifTWijl and Zhao et al. (EEEJ- Unfortunately 
metallicities by Zhao et al. (|200ip could not be consid- 
ered because of an error in their temperature scale, and 
those by McWilliam ( 1990) had to be corrected from a sys- 
tematic trend as described in Paper IV. Complementary 
data were obtained for 88 stars observed with the echelle 
spectrograph ELODIE in February 2003, October 2003, 
and February 2004 at the Observatoire de Haute Provence 
(France) with signal to noise ratios at 550 nm (S/N) rang- 
ing from 150 to 200. We measured and determined their 
radial velocities, [Fe/H] metallicities, and abundances of 
several chemical elements. Metallicities are missing for 7 
remaining stars, among which there are three binaries. 

A detailed description of the atmospheric parameters 
and abundance determination is given in Kovtyukh et al. 
lj2TJU5jl and Mishenina et al. pUTJ5jl . Briefly, the effective 
temperatures were determined with line-depth ratios, a 
technique similar to the one developed by Gray (1994), 
leading to excellent precision of 10-20 K. The surface grav- 
ities logg were determined using two different methods : 
ionisation-balance for iron and fitting the wings of a Ca I 

2 The Hipparcos star positions are expressed in 
the International Celestial Reference System (see 
:/ /www. iers.org/iers/earth/icrs/icrs. html ). 



line profile. For the method of ionisation-balance, we se- 
lected about 100 Fe I and 10 Fe II unblended lines based 
on the synthetic spectra calculations obtained with the 
software STARSP (Tsvmbal. 11555^ . For the profile-fitting 
method, the Ca I line at 6162 A, which is carefully inves- 
tigated in Cayrel et al. 0.996), was used. The gravities 
obtained with these two methods show very good agree- 
ment, as shown in Mishenina et al. (2005). The [Fe/H] 
determination is constrained by the large number of lines 
of neutral iron present in the spectra of giants. The iron 
abundances were determined from the equivalent width of 
lines by applying the program of Kurucz WIDTH9. The 
measurement of the equivalent width of lines was carried 
out with the program DECH20 fGalazutdino\ ll994[) . 

2.2. NGP K giants 

The distant K giant sample was drawn from the Tycho-2 
catalogue (H0g et al., 2000). We applied similar criteria as 
in Paper I to build the list of red clump candidates, just ex- 
tending the limiting apparent magnitudes to fainter stars. 
In summary, we selected stars in two fields close to the 
NGP. The first field (radius 10°, hereafter F10) is centred 
on the NGP, the second one (radius 15°, hereafter F15) is 
centred on the Galactic direction (/ = 35.5°, b = +80°), 
avoiding the Coma open cluster area (a circular field of 
4.5° radius around I = 221°, b = 84°). The total area ef- 
fectively covered by our samples is 720 square degrees. We 
selected stars with 0.9 < B~V < 1.1, in the V magnitude 
range 7.0-10.6 for F10, and 7.0-9.5 for F15 (Johnson mag- 
nitudes transformed from Tycho2 ones) . Known Hipparcos 
dwarfs were rejected. 

A total of 536 spectra were obtained with the 
ELODIE echelle spectrograph at the Observatoire de 
Haute Provence, corresponding to 523 different stars: 347 
in F10 and 176 in F15. The spectra have a median S/N ra- 
tio of 22 at 550 nm. This low S/N is sufficient to estimate 
with good accuracy the stellar parameters, T e ff, gravity, 
and [Fe/H] metallicity, and absolute magnitude My with 
the tgmet method (Katz et al. 11998). as previously de- 
scribed in Paper I. tgmet relies on a comparison by min- 
imum distance of the target spectra to a library of stars 
with well-known parameters, also observed with ELODIE 
(Soubiran et al. Prugniel & Soubiran HUUT^i . Since 

our previous study of clump giants at the NGP, the tg- 
met library has been improved considerably. Many stars 
with well-determined atmospheric parameters and with 
accurate Hipparcos parallaxes have been added to the li- 
brary as reference stars for (T e g, logy, [Fe/H], My). The 
improvment of the tgmet library and the extended sam- 
ple is fully described in Paper IV. Here we just give useful 
characteristics of the extended sample. 

The accuracy of the tgmet results was assessed with 
a bootstrap test on reference stars with very reliable at- 
mospheric parameters and absolute magnitudes. A rms 
scatter of 0.27 was obtained on My and 0.13 on [Fe/H]. 
These values give the typical accuracy of the tgmet re- 
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suits. The scatter on My corresponds to a 13% error in 
distance. As explained in Paper I, the absolute magni- 
tudes from Hipparcos parallaxes of the tgmet reference 
stars are affected by the Lutz-Kelker bias. This causes an 
additional external error that must be taken into account. 
In Paper II, the Lutz-Kelker bias for local clump giants is 
estimated to be —0.09 mag corresponding to a systematic 
overestimation of distances by 4%. We did not attempt 
to correct individual absolute magnitudes of the reference 
stars, because the Lutz-Kelker bias was estimated from 
the luminosity function of the parent population, which is 
unknown for most giants of the TGMET library, except 
for the clump ones. 

In order to test the internal precision of tgmet on 
[Fe/H], we compared the results obtained for the 13 clump 
giants observed twice (Fig.^J. As can be seen, the agree- 
ment is excellent (<r = 0.03 dex). For My, the scatter is 
only 0.08 mag. 

Once the stellar parameters (Toff, logg, [Fe/H], My) 
were determined for each of the 523 target stars of the 
NGP sample, My was used to identify the real red clump 
giants from the dwarfs and subgiants and to compute their 
distances. ELODIE radial velocities, Tycho-2 proper mo- 
tions, and distances were combined to compute 3D ve- 
locities (U, V, W) with respect to the Sun. According to 
typical errors of 0.1 km. s -1 in radial velocities, 15 % in dis- 
tances, 1.4mas.yr _1 in proper motions at a mean distance 
of 470 pc for F10, and 1.2 mas.yr -1 in proper motions at 
a mean distance of 335 pc for F15, the mean errors on the 
two velocity components U and V are 5.6 and 4.0km.s _1 
in F10 et F15, respectively, while it is O.lkm.s -1 on W in 
both fields. 

The 523 target stars are represented in Fig. |2 in 
the plane metallicity - absolute magnitude, with different 
colours for F10 and F15. The F10 field is more contami- 
nated by dwarfs. 



2.3. Binaries 

Binarity can be recognized from the shape of the cor- 
relation function of the spectra with a stellar template. 
None of the 536 target spectra correspond to a clear SB2 
(double-lined spectroscopic binary), but a dozen of them 
present an asymmetric or enlarged profile that could be 
the signature of a companion. Only multi-epoch measure- 
ments could establish the binarity fraction in our sam- 
ple precisely, but it seems to be small as also found by 
Setiawan et al. (2004). This was expected, as for a short 
period binary, the system would merge during the giant 
phase. Moreover, a binary system with a red clump giant 
and a sufficiently bright star, i.e. spectral type earlier than 
K, is expected to have a colour outside our colour selection 
interval. The only star that was really a problem for tg- 
met is TYC 1470-34-1, because of a very enlarged profile, 
which can be due either to a companion or to rotation, so 
this star was removed from the following analysis. 

2.4. Kinematics and metallicity 

We note a deficiency of stars with metallicities [Fe/H] 
< —0.6, although this is not specific to our red clump sam- 
ple (Fig. 0). Our selection criteria does not favour very 
low metallicity stars for which the red clump covers an 
extended colour interval, much larger than our colour in- 
terval selection. This lack of stars with [Fe/H] < —0.60 
exists also in distance-complete samples of dwarfs, for in- 
stance in Reddy et al. I|2003|l . It is also a result obtained 
by Haywood (2001, 2002) when revisiting the metallic- 
ity distribution of nearby G dwarf stars. A more com- 
plete discussion of this aspect of our sample is presented 
in Paper IV. 

For our K z analysis, we rejected stars with 
[Fe/H] metallicities outside the metallicity interval 
[—0.25, +0.10]. TableQlgives the number of stars used for 
the K z determination. 



O. Bienayme et al.: The mass density in the Galactic plane 



5 



Table 1. Number of observed stars and number of iden- 
tified red clump K-giants. 





Hipparcos 


Field 1 


Field2 


Area (sq.deg.) 




309 


410 


V lim 




7-9.5 


7-10.6 


full sample 


203 


176 


347 


red clump My : to 1.3 


203 


124 


204 


Fe cuts -0.25 to +0.10 


152 


67 


100 



Forthy clump stars have [Fe/H] abundances within 
[—0.6, —0.45] and a mean vertical velocity of — 20km.s~ 1 , 
which differs (at 3<r) from the usual or expected 
— 8km. s -1 . The 116 clump stars with [Fe/H] abundances 
within [—0.45, —0.25] have a high velocity dispersion of 
29. 4km. s -1 , and their vertical density distibution de- 
creases slowly over the limited z-extension of our samples. 
Their density decreases too gradually to bring an efficient 
constraint on our present K z determination. Moreover, 
including these stars degrades the analysis by increasing 
the uncertainty on K z drastically. For all these reasons, 
we rejected the lowest metallicity stars from our analysis. 

A small fraction of stars have [Fe/H] > +0.10. These 
stars (24 stars in the local sample and 1 in the NGP cone 
samples) have a relatively low vertical velocity dispersion 
cr^=10km.s _1 ; they also have a correspondingly low scale 
height. In this study, including or rejecting these stars has 
no influence: here they do not constrain the total surface 
mass density of the disk, and have not been included in 
the analysis. 

3. Volume and surface mass density 
determinations 

3.1. The Oort limit 

Thanks to the Hipparcos data, the Galactic potential was 
probed in the first two hundred parsecs from the Sun, 
giving excellent accuracy for determining the Oort limit, 
i.e. the total volume density in the Galactic plane at the 
solar position. 

A first set of studies based on A to F Hipparcos stars 
within a sphere of about 125 pc, have been published 
by Pham (JEHU, Creze et al. (|1998al I1998b|> . and by 
Holmberg & Flynn (2000). Their results all agree within 
1-cr error limits, with differences depending probably on 
the various ways the potential has been parameterized. 
Holmberg & Flynn (2000) model the local potential ac- 
cording to a set of disk stellar populations and also a thin 
disk of gas, while Creze et al. (1998a, 1998b) assume a 
simple quadratic potential. We note that an important 
limitation has been the lack of published radial velocities, 
limiting both the accuracy of the modelling of the kine- 
matics and the possibility of checking the stationary state 
of the various samples used for these analyses indepen- 
dently. From these studies, we will consider (see Holmberg 
& Flynn 2004 and Paper II) that the Oort limit is deter- 



mined as ptotalO = 0) = 0. 100 ± 0.010 M pc- 3 . The Oort 
limit includes the local mass density from both the disk 
and dark halo components. 

Recently Korchagin et al. (2003), using Hipparcos 
data, analyse the vertical potential at slightly larger dis- 
tances. The tracer stars are giant stars that are brighter 
than clump giants, within a vertical cylinder of 200 pc ra- 
dius and an extension of ±400 pc out of the Galactic plane. 
For the dynamical estimate of the local volume density, 
they obtain: p t otai(z = 0) = 0.100 ± 0.005 M pc~ 3 . A 
small improvement could perhaps still be achieved using 
2MASS colour to minimize uncertainties on extinction, dis- 
tances, and vertical velocities. 

3.2. The total disk surface mass density. 
Model 

To analyse the vertical density and velocity distributions 
of our samples and to measure the surface mass density of 
the Galactic disk, we model the /(z, w) distribution func- 
tion, with z the vertical position, w the vertical velocity 
relative to the Local Standard of Rest, and we adjusted 
the free model parameters by a least-square fitting to the 
apparent magnitude star counts a(m) and to the vertical 
velocity dispersions in different magnitude intervals. We 
simultaneously adjusted the model to the data from the 3 
samples. 

The distribution function is modeled as the sum of two 
isothermal components, according to 

/(*,«)= £ exp-(^H^ 2 )^ . (1) 

For the total vertical potential $(z), we used the para- 
metric expression proposed by Kuijken & Gilmore ljl989): 

$(z) ~ S (y/z* + £ 2 - D) + PeS z 2 (2) 

where the potential is related, through the Poisson equa- 
tion, to the vertical distribution of the total volume mass 
density ptotai(^)- The z-integration gives the total surface 
mass density within ± z from the Galactic plane: 

S^kpc = £(< \z\) = , ^ Q2 L 9 +2p cff z. 

Vz + U 

This parametric law models the vertical mass density 
by two density components. One component mimics the 
locally constant density p e g of a round or slightly flattened 
halo. It produces a vertical quadratic potential locally. The 
other component mimics the potential of a flat disk; one 
parameter, Eqi is its surface mass density while the other, 
D, is its half thickness. 

Data 

Compared to our previous work (Paper II), we have in- 
creased the number of observed stars (726 against 387) 
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Table 2. Red clump disk properties. 



[Fe/H] from -0.25 to +0.10 

local relative density velocity dispersion (kms -1 ) 
ci =0.73 ±0.12 (7i = 11.0 ± 1.1 

c 2 = 0.27 ±0.12 cr 2 = 19.7 ±1.7 



and the limiting distances by a factor 1.26. We also ob- 
tained new or first [Fe/H] determinations and built a com- 
plete sample of Hipparcos red clump stars with measured 
metallicities. 

Due to the lack of stars in the range 100-300 pc, our 
samples are not suited to constrain the potential in this 
region efficiently. More generally, pencil beam samples to- 
wards the Galactic poles are not adequate for measuring 
the Oort limit without supplementary assumptions on the 
shape of the vertical potential. Since the Oort limit has 
been previously and accurately measured, we will adopt 
for the volume mass density the value discussed above: 
PtotaiC* = 0) = 0.10 ± 0.01 M 0P c- 3 . 

Fitting the vertical potential with one free parameter: 

A least-square minimization is applied to the observed star 
count distributions and vertical velocity dispersion distri- 
butions (see continuous lines Fig. 13141 and Table We 
fit the data from the three samples, local and distant. We 
fit all the data simultaneously to improve the parameter 
estimation and also the estimate of errors on parameters. 

Individual errors on distances are small for stars from 
the NGP samples, about 13%, but this does not contribute 
to the model uncertainties since the analysed quantities 
(star counts and vertical velocities) are independent of dis- 
tance. Vertical velocities remain slightly affected by uncer- 
tainties on distances through the small contribution of the 
projection of proper motions on the vertical ^-direction. 
The uncertainty on distances, however, reflects the accu- 
racy achieved on absolute magnitude determination for 
distant stars (about 0.27 mag) and our ability to iden- 
tify distant clump stars. Clump star absolute magnitudes 
are normally distributed (see Paper II) around M v = 0.74 
with a dispersion of 0.25, and we select them in the range 
< My < 1.3. Our results also depend on the absolute 
magnitude calibration of nearby giant stars, which is also 
explained in Paper II. 

We may consider that the main source of uncertain- 
ties is the restricted size of the samples. They are split 
into magnitude-bins, and the uncertainty is dominated by 
Poisson fluctuations. For a given bin of magnitude, the er- 
ror on a(m) is y^a(m), and the error on a w is a w /y/2a(m), 
where a(m) is the number of stars in the bin. Errors given 
in Table |3] and elsewhere are deduced from the diagonal 
of the covariance matrix given by the least-square fit. 



density in the Galactic plane 




V 



Fig. 3. Observed red clump counts (circles), respectively 
0.9 and 0.833 magnitude wide bins for the fields 1 and 
2 towards the NGP (local sample and both NGP samples 
are fitted simultaneously) . The continous line indicates the 
best-fit model for the one-parameter potential. 



Results 

We assume that the Galactic dark matter component 
needed to explain the flat rotation curve of our Galaxy 
is spherical and that its density in the solar neighbour- 
hood is p cff = 0.007 M pc~ 3 (Holmberg & Flvnn. |2DUU|| . 
Adopting this value, along with the previously mentioned 
local density /?totai(z = 0), then only one parameter of our 
potential expression is free, since the two parameters Eo 
and D are related through 

ptotei(z = 0) = pefi + S / (2 D) . 

We find D — 260±24pc, and the surface density within 
800pc is found to be E 800 pc = 60 ± 5M pc- 2 . Within 
1.1 kpc, we obtain Ei. ikpc = 64±5 Mqpc -2 . The most re- 
cent determination, £i.ikpc= 71±6 Mqpc -2 , obtained by 
Holmberg & Flynn (2004) is comparable. Their study and 
ours have many similarities, with giant stars in the same 
range of apparent magnitudes, accurate absolute magni- 
tude and distance estimates, similar number of stars and 
apparently similar kinematics, but also many differences: 
NGP versus SGP stars, spectroscopic versus photometric 
distances, clump giants versus a sample with a slightly 
larger range in colours and absolute magnitudes. 

With an estimated mass density of 53 Mqpc -2 for the 
visible matter (Holmberg & Flvnn. I2TJU4) (see Sect. 14.1(1 . 
there is no need a priori to add matter to the model, 
for example, by flattening the dark corona to explain our 
measured local surface density. 
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Fig. 4. Vertical velocity dispersion versus magnitude (see 
Fig.[3J). Velocity dispersions of Hipparcos stars within 
100 pc are plotted at V = 6. 



Table 3. Solutions for the thickness of the vertical poten- 
tial and the total surface mass density. Uncertainties on 
D and £ are ~ 9% 



Pcff 


D 


£800 pc 


Sl.l kpc 


M pc- 3 


pc 


M Q pc^ 


M pc- 2 


0.00 


287 


57 


57 


0.007 


260 


60 


64 


0.01 


249 


61 


67 


0.021 


205 


66 


79 




Fig. 5. The vertical potential (top), the K z force (mid- 
dle), and the total volume density (bottom) for the three 
solutions (p cff = 0,0.01,0.021M Q pc- 3 ) given in Tabled 



The other parameters (relative densities and veloc- 
ity dispersions of stellar components) are explained in 
Table [5] the two fixed parameters are the Sun's position 
above the Galactic plane, zq = 15 pc, and the Sun's vertical 
velocity wq =8 kms -1 . 

Fitting the vertical potential with two free parameters 

Adjusting the model (Eq.|2J with a single free parameter 
gives a satisfying result with small error bars (9 percent), 
but it does not give information on correlations with the 
other (fixed) parameters and does not tell us much about 
the range of possible other solutions that are compatible 
with the observations. But it is known that in practice, 
the change of p e ff does have a significant impact on solu- 
tions: its correlation with others parameters is discussed 
by Gould ITO! . 

Adjusting both Eo and p e ff we find that the best fit is 
obtained with p c ft = 0, while other parameters are very 



close to these given in Table |2 Solutions within 1-cr from 
this best fit give the range 0-0.021 Mqpc -2 for the p e ff 
parameter. 

Table |3 shows some p e g values and solutions that still 
result in acceptable fits to our observations. Thus ac- 
ceptable solutions for Ei.ikpc cover a greater range of 22 
Mqpc -2 than does using a model with just one free pa- 
rameter. This is understandable considering that the best 
determined quantity is the vertical potential (the adjusted 
quantity in Eq. ^) , not the K z force. The derivative of the 
potential gives the surface density, which is proportional 
to the K z force, and similar potentials may produce dif- 
ferent K z forces and surface densities (see Fig. EJ. 

This probably explains the scatter in determinations of 
the surface density in previous works in the 80's and 90's. 
These correlations were already discussed in detail in the 
same context by Gould ( 1990) to explain differences in the 
local surface mass density estimated by various authors. 
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We also note that Eq.|21 parameterizes the vertical po- 
tential and the total vertical mass distribution; however, 
the two apparent components of the r.h.s of Eq.[3 should 
not be strictly read as being one component for the stellar 
and gas disks (with a total surface mass density S ) and 
another one for a round or flattened dark matter halo. 
Only the total vertical mass density distribution (plotted 
at the the bottom of Fig. [SJ is constrained by our data; 
from that figure, we notice that the three plotted solu- 
tions ptotai {z) are quite similar in the range of distances 
to 400 pc. 



For example, in this paper, we have adjusted the verti- 
cal potential with two free parameters; as a consequence, 
the formal errors are larger while the observational 
accuracy of our data is probably better or, at least of 
similar quality, than in the other recent studies. But this 
adjustment with two free parameters allows us to probe 
more realistic and general potentials. After analysing our 
data, we drew the following conclusion: the uncertainty 
on Si.ikpc = 68 ± IIMqpc" 2 (see TableEJ) is still 
large, about sixteen percent, and this must also hold for 
previously published analyses. 



4. Discussion 

An abundant literature and a variety of results exist on the 
vertical K z force and the different methods applied to con- 
strain the total local mass volume density, the Oort limit. 
The oldest papers must be read with care, since system- 
atic bias on distances was more difficult to check before 
the Hipparcos satellite, but the techniques developed and 
comments remain valid. Stellar samples are extracted from 
thin and thick disk populations and are represented as 
the combination of isothermal populations. Kuijken (1991 ) 
proposed an original modeling with a continuum set of 
such isothermal populations. We estimate that the most 
decisive aspect to understand the differences between au- 
thors concerns the various techniques applied in modeling 
the vertical potential. 

For instance, a powerful technique is non-parametric 
modeling of the K z force (Chen et al. 120031 and refer- 
ences in Philip & Lu 1989). These non-parametric K z 
determinations have been achieved without smoothing or 
regularisation (or, for instance, without a positivity con- 
dition on the mass distribution), and the resulting pub- 
lished K z forces show large oscillations that certainly re- 
sult from the small sample sizes. One consequence is that 
these oscillations have no physical interpretation; for in- 
stance, the resulting total mass distribution is not positive 
everywhere. We expect that a conveniently applied regu- 
larisation should be sufficient for making these methods 
more reliable in this context. 

On the other hand, the parametric modeling consists in 
assuming a global shape for the vertical potential. Bahcall 
l|1984|l and recently Holmberg & Flynn H2004j) have as- 
sumed a total disk mass proportional to the observed disks 
of gas and stars, with an extra component of dark matter. 
This last component, proportional to one of the known 
components, is adjusted to constrain the vertical poten- 
tial. Kuijken & Gilmore (1989 ) also proposed a simple an- 
alytical model (the model used in this study) and adjusted 
one of the parameters. 

The advantage of using such a priori knowledge and 
realistic models is to minimize the number of free param- 
eters and to reduce the measured uncertainties. But this 
decrease does not mean that the adjusted parameters have 
been really obtained with better accuracy, since the infor- 
mation on correlations is lost. 



4.1. The visible surface mass density: 

Similar values of £*, the stellar surface mass den- 
sity, at the solar position have been proposed in re- 
cent works: E* ~ 25Mqpc~~ 2 fChabrier. I2001[> . but also 
29M pc- 2 (Holmberg & Flynn, ,2000} and 28M pc" 2 
from the Besangon Galaxy model (Robin, I2003[l . The 
brown dwarf contributions to the surface mass density 
have been estimated: 6Mqpc~ 2 (Holmberg & Flynn, 
l2T)0Hl or 3M pc- 2 (Chabrier, EES)- More difficult is to 
estimate the ISM contribution; Holmberg & Flynn (j^OOO ) 
proposed 13Mqpc~ 2 , but this quantity is uncertain. For 
instance it has been proposed that all the dark matter 
could contribute to the ISM disk component. Adding all 
known contributions, Holmberg & Flynn (2004) propose 
53Mqpc -2 , the value we adopt here. 

4.2. Are our K z solutions compatible with our current 
knowledge of the Galactic rotation curve? 

To answer this question, we simplified and adopted a 
double-exponential density distribution for the Galactic 
disk, including stars and the interstellar medium. We set 
the scale length I — 3kpc, the scale height h = 300 pc, 
and its local surface density £(i?o) = 53Mqpc -2 at the 
solar radius Rq, which includes 40Mqpc -2 for the stel- 
lar contribution and 13Mqpc -2 for the ISM. We ne- 
glected the contributions of the bulge and the stellar halo, 
as these are very small beyond 3kpc from the Galactic 
centre. To maintain a flat Galactic rotation curve, we 
added a Miyamoto spheroid (Miyamoto and Nagai,[l975). 
Adopting Rq = 8.5 kpc and a flat Galactic rotation curve 
V C (R = 5 to 20 kpc) = 220kms _1 , we adjusted the core 
radius (a + b) and the mass of the Miyamoto spheroid to 
find a + b = 9.34 kpc. If the Miyamoto component is spher- 
ical (<z=0kpc), its local density is pd.m. = 0.012 Mqpc -3 . 
Flattening this dark component (a=5kpc, i.e. an axis 
ratio of 0.51), we obtained pd.m. = 0.021 Mqpc~ 3 : the 
exact limit compatible with our K z analysis that still 
gives good fits to Ei.ikpc- Larger flattenings are excluded, 
as, for instance, a= 6.5 kpc gives a 0.34 axis ratio and 
Pd.m. = 0.030 M pc- 3 . 

In conclusion, our data shows that the Galactic dark 
matter can be distributed in a spherical component, but 
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it certainly cannot be distributed in a very flattened disk. 
For instance, if the Galaxy's mass distribution were to- 
tally flat, the local surface density would be as high as 
211M pc- 2 (with V C (R) = 220 km s^ 1 and i? =8.5kpc 
for a Mestel disk, see Binney & TremaineQ.987). In con- 
clusion, there is room to flatten the dark matter halo by a 
maximum factor of about two or three. This agrees with 
the shape of the dark halo that de Boer ( 2005 ) interprets 
from diffuse Galactic EGRET gamma ray excess for ener- 
gies above 1 GeV. 

4.3. Terrestrial impact cratering 

The main topic of this paper is the determination of the 
local surface mass density for that purpose we used 

the most recent determinations of the local volume mass 
density ptotaiO = 0) = 0.10 ± 0.01 M pc~ 3 . Our data 
alone gives a less accurate but close value of ptotai(^ = 
0) = 0.1 ± 0.02 Mqpc -3 . All these recent determinations 
are based on very different samples of stars: A-F dwarfs 
and different types of giants. These samples cover different 
distances, either in a sphere of 125 pc around the Sun, 
inside a cylinder of 800 pc length crossing the Galactic 
plane, or in pencil beams up to 1.1 kpc from the plane. 
All these determinations converge towards the same value 
of the local volume mass density, 0.10 Mqpc~ 3 , implying 
that the half period of the vertical oscillations of the Sun 
through the Galactic plane is 42±2Myr. 

A 26Myr periodicity in epochs of major mass ext- 
inction of species was found by Raup & Sepkoski (|1986|l . 
and a cycle of 28.4 Myr in the ages of terrestrial impact 
craters was found by Alavarez & Muller l|1984fl . The peri- 
odicity in these catastrophes is disputed, however (see for 
instance Jetsu & Pelt. 12000)1 . The spectral analysis of the 
periodicity hypothesis in cratering records shows, in the 
most recently published works, possible or significant pe- 
riods: 33 ±4 Myr (Rampino & Stothers, 1984), 33 ± lMyr 
(Stothers. ll998[) . and more recently 16.1 Myr and 34.7Myr 
by Moon et al. (|2003|) . We note that some authors estimate 
that periodicities could result from a spurious "human- 
signal" such as rounding (Jetsu & Pelt. 120077)) . Recently, 
Yabushita l|2002l l2TTO)l claims a periodicity of 37.5 Myr 
and considers that the probability of deriving this period 
by chance on the null hypothesis of a random distribu- 
tion of crater ages is smaller than 0.10. From the width of 
the peaks in their periodograms, we estimate their period 
accuracy to be ±2 Myr. 

It has been frequently proposed that the period 
of high-impact terrestrial cratering would be directly 
linked to the crossing of the solar system through the 
Galactic plane where the giant molecular clouds are 
concentrated. However, the above-mentioned periods 
in the range 33-37 Myr would correspond to the half 
period of oscillation of the Sun only if the total local 
mass density were 0.15 Mqpc -3 , a high value of the local 
density measured in the 1980's by a few authors. Recent 
and accurate ptotal(z — 0) measurements now imply a 



42±2 Myr half period of oscillation of the Sun that can 
no longer be related to possible periods of large impact 
craters and mass extinction events. 



5. Summary 

• Adopting the previous determination of the total lo- 
cal volume density ptotei(z = 0) = 0.10 Mqpc~ 3 in the 
Galactic disk, we modelled the vertical disk potential and 
mass distribution through the constraints of a sample of 
red clump stars towards the NGP with measured dis- 
tances, velocities, and [Fe/H] abundances. 

• Our simplest model, including a spherical dark corona, 
shows that there is no need for extra dark matter in the 
Galactic disk to explain the vertical gravitational poten- 
tial. The total surface mass density, at the Solar position, 
is found to be Si. i kpc = 64 ± 5 Mqpc -2 , compared to 
53 Mqpc -2 for the visible matter (see section l4~Tl . which 
is sufficient to explain both the observed amount of visible 
matter and the local contribution of a round dark matter 
halo (15 M pc- 2 ). 

• With a two-parameter model and some flattening of 
the dark matter halo, we obtained a larger extent of ac- 
ceptable solutions: Ei.ikpc = 57-79 Mqpc -2 , £o.8k pc = 57- 
66Mqpc~ 2 . Flattening larger than about 2 or 3 is ex- 
cluded by analysis of our red clump giants. 

We note that a flattening by a factor two of the dark 
corona could contradict the other recent constraint ob- 
tained by Ibata et al. ( 2001), which is based on the non- 
precessing orbit of the Sagittarius stream. However, while 
their result concerns the outer halo, our analysis is only 
sensitive to the inner one. 

On the other hand, it is not possible to flatten the 
dark matter halo without increasing its own contribution 
to the local volume mass density too much. The recent dy- 
namical estimates based on Hipparcos data of the total 
volume mass density gives 0.10 Mqpc -3 . This includes the 
known stellar local mass density, 0.046 Mqpc -3 , and the 
gas volume mass density, 0.04 Mqpc -3 (see Palasi 1998 
Holmberg & Flynn IM7U1 Chabrier I27M1 I27JU1 12HU1 or 
Kaberla 2003). This leaves room for only 0.014 Mqpc 
for the dark matter. This also implies that the halo cannot 
be flattened more than a factor ^two, unless the volume 
mass density of the gas, the weakest point in the K z anal- 
ysis, has been strongly overestimated. 

• As a by-product of this study we determined the half 
period of oscillation of the Sun through the Galactic plane, 
42 ± 2 Myr, which cannot be related to the possible period 
of large terrestrial impact craters ~ 33-37 Myr. 
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